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Abstract 

The physical implementation of the quantum Control-Not gate for a two-spin sys- 
tem is investigated numerically. The concept of a generalized quantum Control-Not 
gate, with arbitrary phase shift, is introduced. It is shown that a resonant 7r-pulse 
provides a simple example of a generalized quantum Control-Not gate. 
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Recently, the field of quantum computation has experienced remarkable progress. 
(See, for example, the review M). Important achievements include a practical imple- 
mentation of quantum computing 0, a quantum algorithm for prime factorization ||, 
a quantum error corrections code M and an algorithm for pattern recognition ||. It 
was shown in || that two qubits quantum Control-Not (CN) gates in combination 
with one qubit rotations can provide all the quantum logic gate. Because of this, sig- 
nificant attention has been directed toward the implementation of quantum CN gates 
[0]"IH3- The simplest practical implementation of quantum CN logic gates is based 



on the application of a 7r-pulse to a chain of weakly interacting two-level atoms (ions, 
spins) which have different frequencies of transition ||. A qubit can be spanned by 
the ground and long-lived exited state of an atom. Weak interaction between atoms 
provides the implementation of the quantum CN gate. 

In our previous paper [jTCfl , we explored this idea for the system of two nuclear 
spins connected by the Ising interaction. We investigated the evolution of the complex 
amplitudes, Cjfc, of the wave function, 

2 

^(*) = E *k(t)\ik >, (1) 

i,k=0 

under the action of a resonant 7r-pulse. The first number % in \ik > indicates the state 
of the control spin (qubit); the second number k indicates the state of the target spin 
(qubit). Our numerical calculations successfully confirmed the idea |J concerning the 
behavior of the moduli of the quantum complex amplitudes. For quantum computa- 
tions, the dynamics of phases of complex amplitudes Cik is also important. Indeed, 
quantum interference effects depends on both, moduli and phases. 

In this paper, we investigate the dynamics of the quantum CN gate taking into 
consideration the behavior of moduli and phases of the complex amplitudes, c^. We 
consider the system with the Hamiltonian |1C| . 



(2) 



This Hamiltonian describes a four level quantum system with the frequencies: uj n ± J. 
The idea of implementation of a quantum CN gate is the following. Assume that 
one applies to the system (2) a 7r-pulse with the frequency, to = 0J2 — J ■ This pulse 
drives the second (target) spin from the ground state |0 > to the exited state |1 > (or 
vice versa) if the first (control) spin is in exited state |1 >. We recall that the quantum 
CN gate is commonly defined as, 

CN = |00 >< 00| + |01 >< 01 1 + |10 >< 1 1 1 + |11 >< 10|. (3) 

It follows from (3) that any superposition state, 

V = c 00 (0)|00 > +c i(0)|01 > +cio(0)|10 > +c n (0)|ll >, (4) 
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after the action of the CN operator (3), is transformed to the state, 

^ = c 00 (0)|00 > +coi(0)|01 > +ci (0)|ll > +cn(0)|10 > . (5) 
We now introduce the generalized quantum CN gate (GCN gate), 

GCN(A<f ik ) = exp(iAv%))|00 >< 00| + exp(iAy9 i)|01 >< 01|+ (6) 

exp(iA<pn)|10 >< 1 1 1 + exp(iA^i )|ll >< 10|. 

Similar to a "pure" CN gate (3), the GCN gate changes the state of the target qubit 
if the control qubit is in the state |1 >. In addition, it changes the phases of the 
complex amplitudes c^. (Up to an insignificant common phase factor, we can put one 
of the phases Aip ik to be zero.) For quantum computations, it is important to know 
the values of Atf ik , otherwise one can not describe properly the action of a GCN gate. 

Consider the action of a 7r-pulse with frequency w on a system with the Hamilto- 
nian (2). We add to the Hamiltonian (2) the interaction with a circularly polarized 
transverse magnetic field, 

Him = --X! a « (exp(-^)<7~ + exp(?cjt)(T+) , (7) 

71=1 

and substitute the wave function (1) into the Schrodinger equation. To obtain the 
equations of motion with constant coefficients, we use substitution which is equivalent 
to a transition to a rotating frame: 

coo -> c 00 exp[i(ujt + tp)], (8) 

c i -> coi exp(iy?) 

Cio -> c 10 exp(i(p), 

en — > en exp[i(—ut + if)]. 

In (8), tp — (f{t) is a common phase (which can be chosen arbitrary, to simplify the 
equations). We shall choose, 

uj = oj 2 — J, fit) = (cj 2 — J — Ui)t/2. (9) 

Then, we derive the following equations (with time-independent coefficients) for the 
amplitudes c^, 

-2ic 00 = -2{u 2 -ui- 2J)c 00 + aiCio + a 2 c 01 , (10) 
-2ic 01 = -2(cj 2 - ^i)c i + axCu + a 2 c 00 , 
-2ic w = aicoo + a 2 cu, 
—2it\\ = aicoi + a 2 c w . 
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In matrix representation Eqs (10) have the form, 

2icj^ Bikj m Cj m ^ (H) 

where B^ m is a matrix with time-independent elements. The last two equations in 
(10) describe the resonant transition between the state 1 1 > and |11 >. The first two 
equations in (10) describe the nonresonant dynamics of the lower energy states, 1 00 > 
and 1 1 >. 

For the values of parameters, 

uj\ = 5ui2, ^2 = 100, J = 5, oj = u>2 — J = 95, a\ = 0.5, a>2 = 0.1, (12) 

we obtain the dependence of the complex amplitudes on time during the action of the n- 
pulse. (The characteristic dimensional parameters can be chosen: ui2/2ir = 100 MHz, 
Ldi/2-K = 500 MHz.) Note that the effective angle of rotation of the second spin 
under the action of the electromagnetic pulse is slightly larger than g^t, where r is 
the duration of a 7r-pulse. The reason is a weak indirect excitation of the resonant 
transition |10 >^ |11 > via the nonresonant (first) spin. The terms aiCoo and aiCoi in 
the last two equations in (10) are responsible for this effect. 

The free evolution of the two-spin system in the rotating frame (ai = a 2 = in 
equations (10)) is described by the amplitudes, 

coo(t) = c 00 (0) exp[-i(^ 2 - o>i - 2J)t], (13) 

Coi(t) = c i(0) exp[-i(w 2 - u^t], 

cio(i) = Cio(0), c u (t) = cn(0). 

To eliminate the phase factor, corresponding to the free evolution, we will discuss below 
the dynamics of the coefficients c' ik , 

c do(0 = c ooW exp[i(u 2 - u x - 2J)f], (14) 

c diW = c oiW exp[?(cj 2 - 

c 'io(0 = °io(t), c' u (t) = c n (t). 

In Fig. la the time dependence of the real part of c' u , Rec' u (t), and the imaginary 
part of c' 10 (t), Imc' 10 (t), are shown for the initial conditions, 

4(0) = !, c; m (0) = 0, (j,m)^(l,l). (15) 

One can see the monotonic decrease of i?ec' 11 (t) and increase of Imc' 10 (t). At the end 
of the 7r-pulse, c' 10 (t) = ic' u (0). The values of Imc' u (t) and Rec' 10 (t) are negligible as 
well as the values |cq (£)| and |c 01 (£)|. This evolution describes the transformation, 

|ll>->i|10>. (16) 
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In Fig. lb the analogous dependences are shown for the initial conditions 

4,(0) = 1, c; m (0) = 0, (j,m)^(l,0). (17) 

At the end of the 7r-pulse, we have: c' 11 (r) = ic' 10 (0) which corresponds to the trans- 
formation, 

|10>->i|ll>. (18) 

For the initial conditions, 

4(0) = 1, c; m (0) = 0, (j,m)^(0,l), (19) 

(which corresponds to the population of nonresonant level 1 01 >), the amplitudes c' ik {t) 
do not change within an accuracy of 1CT 3 . The same is true for the initial conditions, 

c' 00 (0) = 1, c ; m (0) = (0,0), (j,m)^(0,0). (20) 

In Fig. 2 the time dependence of the amplitudes is shown for the superpositional initial 

state, 

4(0) = v/3/10, 4(0) = l/>/5, ci (0) = 1/A c' n (0) = l/v^G. (21) 

One can see that at the end of the 7r-pulse, the amplitudes take the following values, 

c' n (r)=< (0), 4(r)=<(0). (22) 

(The nonresonant amplitudes c' 0l (t) and c' 00 (t) do not change within an accuracy of 
10~ 3 ). One can conclude that a 7r-pulse with the frequency uj 2 — J implements a GCN 

gate, 

CN(0, 0, tt/2, tt/2) = |00 >< 00 1 + |01 >< 01] +i|10 >< 11| +i\U >< 10|. (23) 

GCN gate (23) can not be used for implementation of a "pure" CN gate (3). 

We also checked an opportunity for implementation of "pure" quantum CN gate 
for the amplitudes Cik{t) in (10), which include fast oscillations of the free evolution. 
This is, probably, not important for quantum computations, but may be interesting 
for experimental investigation of nonresonant quantum states in this system. We have 
found that by changing only the amplitude of a 7r-pulse, we can not get a "pure" CN 
gate. It was necessary to change slightly one of the parameters {jJ\-,oj 2 or J in (12) to 
get a "pure" CN gate. 

In Fig. 3 the action of the "pure" CN gate is demonstrated for the values of 
parameters, 

Ul = 500.06, cu 2 = 100, J = 5, a 2 = 0.10016, ai = (24) 
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and for "digital" initial conditions, 

^(0) = 1, c jm (0) = 0, (j,m) ^ (i,k), (25) 

where the indices ik take the values 00,01, 10 and 11 from the top to the bottom of 
Fig. 3. Because of the fast oscillations of nonresonant amplitudes, we show only the 
evolution of the amplitudes only near the end of the 7r-pulse. One can see that for any 
initial "digital" condition, the only nonzero amplitude at the end of a 7r-pulse has the 
value "i" . In Fig. 4 the action of the same gate is shown for the superpositional initial 
conditions (21). One can see that at the end of the 7r-pulse, 

coo(t) = icoo(0), coi(r) =icoi(0), c w (r) = ic u (0), c u (t) = ic 10 (0) . (26) 

Thus, up to insignificant common phase factor, tt/2, one has the "pure" CN gate (3). 

We are grateful to D.K. Ferry for fruitful discussions. This research was supported 
in part by the Linkage Grant 93-1602 from the NATO Special Programme Panel on 
Nanotechnology. Work at Los Alamos was supported by the Defense Advanced Re- 
search Projects Agency. 
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Figure Captions 



Fig. 1. Time evolution of the amplitudes c' ik under the action of a 7r-pulse, for the 
initial conditions (a) (15) and (b) (17). In (a) curve (1) corresponds to Rec' u (t) and 
curve (2) corresponds to Imc' 10 (t). In (b) curve (1) corresponds to Rec' l0 (t) and curve 
(2) corresponds to Imc' n (t). The vertical arrows show the beginning and the end of 
the 7r-pulse. 

Fig. 2. Time evolution of the amplitudes c' ik under the action of a 7r-pulse, for the 
superpositional initial conditions (21). In (a) curve (1) corresponds to Rec' w (t) and 
curve (2) corresponds to Imc' 10 (t). In (b) curve (1) corresponds to Rec' u (t) and curve 
(2) corresponds to Jmc' 11 (t). The vertical arrows show the beginning and the end of 
the 7r-pulse. 

Fig. 3 Imaginary parts of the amplitudes Qfc(t) near the end of a 7r-pulse for dig- 
ital initial conditions (25), and the values of parameters (24). 

Fig. 4. Imaginary parts of the amplitudes Qfc(t) near the end of a 7r-pulse for the 
superpositional initial conditions (21), and the values of parameters (24). 
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